The Power of SOFIA/FORCAST in Estimating Internal Luminosities of Low Mass Class 0/I Protostars

With the Stratospheric Observatory for Infrared Astronomy (SOFIA) routinely operating science flights, we demonstrate that observations with the Faint Object infraRed CAmera for the SOFIA Telescope (FORCAST) can provide reliable estimates of the internal luminosities, $L_{\rm int}$, of protostars. We have developed a technique to estimate $L_{\rm int}$ using a pair of FORCAST filters: one "short-wavelength" filter centered within 19.7-25.3 $\mu$m, and one "long-wavelength" filter within 31.5-37.1 $\mu$m. These $L_{\rm int}$ estimates are reliable to within 30-40% for 67% of protostars and to within a factor of 2.3-2.6 for 99% of protostars. The filter pair comprised of F25.3$\mu$m and F37.1$\mu$m achieves the best sensitivity and most constrained results. We evaluate several assumptions that could lead to systematic uncertainties. The OH5 dust opacity matches observational constraints for protostellar environments best, though not perfectly; we find that any improved dust model will have a small impact of 5-10% on the $L_{\rm int}$ estimates. For protostellar envelopes, the TSC84 model yields masses that are twice those of the Ulrich model, but we conclude this mass difference does not significantly impact results at the mid-infrared wavelengths probed by FORCAST. Thus, FORCAST is a powerful instrument for luminosity studies targeting newly discovered protostars or suspected protostars lacking detections longward of 24 $\mu$m. Furthermore, with its dynamic range and greater angular resolution, FORCAST may be used to characterize protostars that were either saturated or merged with other sources in previous surveys using the Spitzer Space Telescope or Herschel Space Observatory.Comment: 17 pages, 9 figures. Accepted for publication in Ap

The Spitzer c2d Survey of Nearby Dense Cores. IX. Discovery of a Very Low Luminosity Object Driving a Molecular Outflow in the Dense Core L673-7

We present new infrared, submillimeter, and millimeter observations of the dense core L673-7 and report the discovery of a low-luminosity, embedded Class 0 protostar driving a molecular outflow. L673-7 is seen in absorption against the mid-infrared background in 5.8, 8, and 24 micron Spitzer images, allowing for a derivation of the column density profile and total enclosed mass of L673-7, independent of dust temperature assumptions. Estimates of the core mass from these absorption profiles range from 0.2-4.5 solar masses. Millimeter continuum emission indicates a mass of about 2 solar masses, both from a direct calculation assuming isothermal dust and from dust radiative transfer models constrained by the millimeter observations. We use dust radiative transfer models to constrain the internal luminosity of L673-7, defined to be the luminosity of the central source and excluding the luminosity from external heating, to be 0.01-0.045 solar luminosities, with 0.04 solar luminosities the most likely value. L673-7 is thus classified as a very low luminosity object (VeLLO), and is among the lowest luminosity VeLLOs yet studied. We calculate the kinematic and dynamic properties of the molecular outflow in the standard manner, and we show that the expected accretion luminosity based on these outflow properties is greater than or equal to 0.36 solar luminosities. The discrepancy between this expected accretion luminosity and the internal luminosity derived from dust radiative transfer models indicates that the current accretion rate is much lower than the average rate over the lifetime of the outflow. Although the protostar embedded within L673-7 is consistent with currently being substellar, it is unlikely to remain as such given the substantial mass reservoir remaining in the core.Comment: 19 pages, 14 figures. Accepted by Ap

The Dynamical State of Barnard 68: A Thermally Supported, Pulsating Dark Cloud

We report sensitive, high resolution molecular-line observations of the dark cloud Barnard 68 obtained with the IRAM 30-m telescope. We analyze spectral-line observations of C18O, CS(2--1), C34S(2--1), and N2H+(1--0) in order to investigate the kinematics and dynamical state of the cloud. We find extremely narrow linewidths in the central regions of the cloud. These narrow lines are consistent with thermally broadened profiles for the measured gas temperature of 10.5 K. We determine the thermal pressure to be a factor 4 -- 5 times greater than the non-thermal (turbulent) pressure in the central regions of the cloud, indicating that thermal pressure is the primary source of support against gravity in this cloud. This confirms the inference of a thermally supported cloud drawn previously from deep infrared extinction measurements. The rotational kinetic energy is found to be only a few percent of the gravitational potential energy, indicating that the contribution of rotation to the overall stability of the cloud is insignificant. Finally, our observations show that CS line is optically thick and self-reversed across nearly the entire projected surface of the cloud. The shapes of the self-reversed profiles are asymmetric and are found to vary across the cloud in such a manner that the presence of both inward and outward motions are observed within the cloud. Moreover, these motions appear to be globally organized in a clear and systematic alternating spatial pattern which is suggestive of a small amplitude, non-radial oscillation or pulsation of the outer layers of the cloud about an equilibrium configuration.Comment: To appear in the Astrophysical Journal; 23 pages, 8 figures; Manuscript and higher resolution images can be obtained at http://cfa-www.harvard.edu/~ebergin/pubs_html/b68_vel.htm

The Spitzer c2d Survey Of Nearby Dense Cores. X. Star Formation In L673 And Cb188

L673 and CB188 are two low-mass clouds isolated from large star-forming regions that were observed as part of the Spitzer Legacy Project "From Molecular Clouds to Planet Forming disks" (c2d). We identified and characterized all the young stellar objects (YSOs) of these two regions and modeled their spectral energy distributions (SEDs) to examine whether their physical properties are consistent with values predicted from the theoretical models and with the YSO properties in the c2d survey of larger clouds. Overall, 30 YSO candidates were identified by the c2d photometric criteria, 27 in L673 and 3 in CB188. We confirm the YSO nature of 29 of them and remove a false Class III candidate in L673. We further present the discovery of two new YSO candidates, one Class 0 and another possible Class I candidate in L673, therefore bringing the total number of YSO candidates to 31. Multiple sites of star formation are present within L673, closely resembling other well-studied c2d clouds containing small groups such as B59 and L1251B, whereas CB188 seems to consist of only one isolated globule-like core. We measure a star formation efficiency (SFE) of 4.6%, which resembles the SFE of the larger c2d clouds. From the SED modeling of our YSO sample we obtain envelope masses for Class I and Flat spectrum sources of 0.01-1.0 M-circle dot. The majority of Class II YSOs show disk accretion rates from 3.3 x 10(-10) to 3 x 10(-8) M-circle dot yr(-1) and disk masses that peak at 10(-4) to 10(-3) M-circle dot. Finally, we examined the possibility of thermal fragmentation in L673 as the main star-forming process. We find that the mean density of the regions where significant YSO clustering occurs is of the order of similar to 10(5) cm(-3) using 850 mu m observations and measure a Jeans Length that is greater than the near-neighbor YSO separations by approximately a factor of 3-4. We therefore suggest that other processes, such as turbulence and shock waves, may have had a significant effect on the cloud's filamentary structure and YSO clustering.University of SouthamptonNASA 1279198, 1288806, 1365763Jet Propulsion Laboratory, California Institute of TechnologyAstronom

Observational Constraints on Submillimeter Dust Opacity

Infrared extinction maps and submillimeter dust continuum maps are powerful probes of the density structure in the envelope of star-forming cores. We make a direct comparison between infrared and submillimeter dust continuum observations of the low-mass Class 0 core, B335, to constrain the ratio of submillimeter to infrared opacity (κ_(smm)/κ_(ir)) and the submillimeter opacity power-law index (κ ∝ λ–β). Using the average value of theoretical dust opacity models at 2.2 μm, we constrain the dust opacity at 850 and 450 μm. Using new dust continuum models based upon the broken power-law density structure derived from interferometric observations of B335 and the infall model derived from molecular line observations of B335, we find that the opacity ratios are ^κ_(850)_κ_(2.2) = (3.21 - 4.80)^(+0.44)_(-0.30) x 10^(-4) ^κ_(450)_κ(2.0) = (12.8-24.8)^(+2.4)_(-1.3) x 10^(-4) with a submillimeter opacity power-law index of β_(smm) = (2.18-2.58)^(+0.30)_(–0.30). The range of quoted values is determined from the uncertainty in the physical model for B335. For an average 2.2 μm opacity of 3800 ± 700 cm^2 g^(–1), we find a dust opacity at 850 and 450 μm of κ_(850) = (1.18-1.77)^9+0.36)_(–0.24) and κ_(450) = (4.72-9.13)^(+1.9)_(–0.98) cm^2 g^(–1) of dust. These opacities are from (0.65-0.97)κ^(OH5)_(850) of the widely used theoretical opacities of Ossenkopf and Henning for coagulated ice grains with thin mantles at 850 μm

The Spitzer c2d Survey Of Nearby Dense Cores. XI. Infrared And Submillimeter Observations Of CB130

We present new observations of the CB130 region composed of three separate cores. Using the Spitzer Space Telescope, we detected a Class 0 and a Class II object in one of these, CB130-1. The observed photometric data from Spitzer and ground-based telescopes are used to establish the physical parameters of the Class 0 object. Spectral energy distribution fitting with a radiative transfer model shows that the luminosity of the Class 0 object is 0.14-0.16 L-circle dot, which is low for a protostellar object. In order to constrain the chemical characteristics of the core having the low-luminosity object, we compare our molecular line observations to models of lines including abundance variations. We tested both ad hoc step function abundance models and a series of self-consistent chemical evolution models. In the chemical evolution models, we consider a continuous accretion model and an episodic accretion model to explore how variable luminosity affects the chemistry. The step function abundance models can match observed lines reasonably well. The best-fitting chemical evolution model requires episodic accretion and the formation of CO2 ice from CO ice during the low-luminosity periods. This process removes C from the gas phase, providing a much improved fit to the observed gas-phase molecular lines and the CO2 ice absorption feature. Based on the chemical model result, the low luminosity of CB130-1 is explained better as a quiescent stage between episodic accretion bursts rather than being at the first hydrostatic core stage.NASA 1224608, 1288664, 1407, NNX07AJ72G, 1279198, 1288806, 1342425NSF AST-0607793, AST-0708158Korea government (MEST) 2009-0062866Ministry of Education, Science and Technology 2010-0008704Astronom

Comparing Star Formation on Large Scales in the c2d Legacy Clouds: Bolocam 1.1 mm Dust Continuum Surveys of Serpens, Perseus, and Ophiuchus

We have undertaken an unprecedentedly large 1.1 millimeter continuum survey of three nearby star forming clouds using Bolocam at the Caltech Submillimeter Observatory. We mapped the largest areas in each cloud at millimeter or submillimeter wavelengths to date: 7.5 sq. deg in Perseus (Paper I), 10.8 sq. deg in Ophiuchus (Paper II), and 1.5 sq. deg in Serpens with a resolution of 31", detecting 122, 44, and 35 cores, respectively. Here we report on results of the Serpens survey and compare the three clouds. Average measured angular core sizes and their dependence on resolution suggest that many of the observed sources are consistent with power-law density profiles. Tests of the effects of cloud distance reveal that linear resolution strongly affects measured source sizes and densities, but not the shape of the mass distribution. Core mass distribution slopes in Perseus and Ophiuchus (alpha=2.1+/-0.1 and alpha=2.1+/-0.3) are consistent with recent measurements of the stellar IMF, whereas the Serpens distribution is flatter (alpha=1.6+/-0.2). We also compare the relative mass distribution shapes to predictions from turbulent fragmentation simulations. Dense cores constitute less than 10% of the total cloud mass in all three clouds, consistent with other measurements of low star-formation efficiencies. Furthermore, most cores are found at high column densities; more than 75% of 1.1 mm cores are associated with Av>8 mag in Perseus, 15 mag in Serpens, and 20-23 mag in Ophiuchus.Comment: 32 pages, including 18 figures, accepted for publication in Ap

The Spitzer c2d Survey of Nearby Dense Cores. V. Discovery of a VeLLO in the "Starless" Dense Core L328

This paper reports the discovery of a Very Low Luminosity Object (VeLLO) in the "starless" dense core L328, using the Spitzer Space Telescope and ground based observations from near-infrared to millimeter wavelengths. The Spitzer 8 micron image indicates that L328 consists of three subcores of which the smallest one may harbor a source, L328-IRS while two other subcores remain starless. L328-IRS is a Class 0 protostar according to its bolometric temperature (44 K) and the high fraction ~72 % of its luminosity emitted at sub-millimeter wavelengths. Its inferred "internal luminosity" (0.04 - 0.06 Lsun) using a radiative transfer model under the most plausible assumption of its distance as 200 pc is much fainter than for a typical protostar, and even fainter than other VeLLOs studied previously. Note, however, that its inferred luminosity may be uncertain by a factor of 2-3 if we consider two extreme values of the distance of L328-IRS (125 or 310 pc). Low angular resolution observations of CO do not show any clear evidence of a molecular outflow activity. But broad line widths toward L328, and Spitzer and near-infrared images showing nebulosity possibly tracing an outflow cavity, strongly suggest the existence of outflow activity. Provided that an envelope of at most ~0.1 Msunis the only mass accretion reservoir for L328-IRS, and the star formation efficiency is close to the canonical value ~30%, L328-IRS has not yet accreted more than 0.05 Msun. At the assumed distance of 200 pc, L328-IRS is destined to be a brown dwarf.Comment: 29 pages, 8 figures, 1 table, to be published in Astrophysical Journa

Bolocam Survey for 1.1 mm Dust Continuum Emission in the c2d Legacy Clouds. II. Ophiuchus

We present a large-scale millimeter continuum map of the Ophiuchus molecular cloud. Nearly 11 square degrees, including all of the area in the cloud with visual extinction more than 3 magnitudes, was mapped at 1.1 mm with Bolocam on the Caltech Submillimeter Observatory (CSO). By design, the map also covers the region mapped in the infrared with the Spitzer Space Telescope. We detect 44 definite sources, and a few likely sources are also seen along a filament in the eastern streamer. The map indicates that dense cores in Ophiuchus are very clustered and often found in filaments within the cloud. Most sources are round, as measured at the half power point, but elongated when measured at lower contour levels, suggesting spherical sources lying within filaments. The masses, for an assumed dust temperature of 10 K, range from 0.24 to 3.9 solar masses, with a mean value of 0.96 solar masses. The total mass in distinct cores is 42 solar masses, 0.5 to 2% of the total cloud mass, and the total mass above 4 sigma is about 80 solar masses. The mean densities in the cores are quite high, with an average of 1.6 x 10^6 per cc, suggesting short free-fall times. The core mass distribution can be fitted with a power law with slope of 2.1 plus or minus 0.3 for M>0.5 solar masses, similar to that found in other regions, but slightly shallower than that of some determinations of the local IMF. In agreement with previous studies, our survey shows that dense cores account for a very small fraction of the cloud volume and total mass. They are nearly all confined to regions with visual extinction at least 9 mag, a lower threshold than found previously.Comment: 47 pages, 16 figures, accepted for Ap